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STRESS  INTENSITY  FACTORS  FOR  CRACKS  AT  LOADED  BOLES 

by 

D.  P.  Rooke 
J.  Tweed* 


SUMMARY 


*  A  method  is  described  for  obtaining  the  opening-mode  and  sliding-mode  stress 
intensity  factors  for  the  tip  of  a  crack  at  the  edge  of  a  circular  hole  with 
arbitrary  loads  on  its  perimeter.  The  method  involves  the  solutions  of  a  singular 
integral  equation  which  arises  from  expressing  the  equations  of  elasticity  in 
terms  of  Mellin  transforms.  In  developing  the  method,  four  particular  load  dis¬ 
tributions  are  considered:  a  uniform  tensile  stress  remote  from  the  holej  a  point 
load  on  the  perimeter  of  the  hole;  uniform  pressure  on  an  arc  of  the  perimeter; 
a  cosine  distribution  of  pressure  on  the  perimeter.  It  is  shown  that  the  stress 
intensity  factor  for  a  short  crack  is  strongly  dependent  on  the  distribution  of 
load  around  the  hole;  this  has  important  implications  for  the  estimation  of 
fatigue  lifetimes.  A  procedure  is  indicated  for  using  some  of  the  results  to 
study  arbitrary  loads  on  the  hole  perimeter. — 
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t  INTRODUCTION 

Cracks  in  structural  components  can  grow  under  in-service  loads  and  lead  to  failure 
of  the  component.  In  order  to  ensure  safety  and  reliability  it  is  necessary  to  know  both 
the  residual  strength  of  a  cracked  structure  and  the  rates  at  which  cracks  grow  in 
fatigue;  both  depend  upon  the  stress  field  at  the  crack  tip,  which  is  characterised  by 
the  stress  intensity  factor. 

Because  holes  are  stress  concentrators,  cracks  frequently  start  at  the  edges  of 
holes  in  the  structure.  Rooke  and  Tweed'  have  evaluated  the  opening  mode  stress  intens¬ 
ity  factors  Kj  for  cracks  at  the  edges  of  circular  holes  in  sheets  subjected  to  remote 
tensile  stresses  which  result  in  only  normal  stresses  along  the  crack  line;  the  sliding 
mode  stress  intensity  factors  were  therefore  zero.  However,  in  many  practical 

components  the  service  loads  produce  shear  stresses  along  the  crack  line,  and  a  more 
general  method  of  analysis  is  required  to  evaluate  the  stress  intensity  factors  (K^j 

will  now  be  non-zero,  but  generally  small  compared  with  K.  ).  Such  a  method  is  now 

.  2  3  1 

available  *  ;  it  can  be  used  to  obtain  both  and  for  a  crack  at  the  edge  of  a 

circular  hole  in  a  sheet  subjected  to  arbitrary  loading. 

In  this  Report  the  aim  was  to  model  two  cosmon  practical  configurations,  namely  a 
crack  at  the  edge  of  a  hole  in  a  sheet  subjected  to  arbitrary  loadings  remote  from  the 
crack,  and  a  crack  at  a  hole  subjected  to  arbitrary  loadings  on  the  perimeter  of  the 
hole.  Four  particular  load  distributions  are  considered,  see  Fig  1:  they  are 

(a)  a  uniform  uniaxial  tensile  stress  remote  from  the  crack  acting  in  a  direction  at  an 
arbitrary  angle  to  the  crack  line; 

(b)  a  localised  radial  force  at  the  perimeter  of  the  hole  acting  at  an  arbitrary  angle 
to  the  crack  line; 

(c)  a  uniform  pressure  acting  over  an  arc  of  the  hole  perimeter  such  that  the  resultant 
force  on  the  edge  of  the  hole  acts  in  a  direction  at  right  angles  to  the  crack 
line; 

(d)  a  variable  pressure  (cosine  distribution)  acting  over  one  half  of  the  hole  peri¬ 
meter  such  that  the  resultant  force  on  the  edge  of  the  hole  acts  in  a  direction  at 
right-angles  to  the  crack  line. 

Stress  intensity  factors  for  cracks  in  structures  subjected  to  loadings  in  more 
than  one  direction  can  be  obtained  from  the  results  for  case  (a)  by  using  the  principle 
of  superposition.  The  stress  intensity  factor  for  the  single  radial  point  fo.ce  in 
case  (b)  can  be  used  as  a  Green's  function  to  obtain  results  for  any  arbitrary  pressure 
acting  on  the  hole  perimeter.  Two  particular  pressure  distributions,  cases  (c)  and  (d), 
have  been  studied  in  detail,  since  these  distributions  have  often  been  assumed  in 
attempts  to  model  the  load  transfer  in  both  pin-loaded  lugs  and  fastener  holes. 

The  results  obtained  in  this  Report  show  that  the  stress  intensity  factors,  Kj 
and  ,  for  short  cracks  depend  significantly  upon  the  distribution  of  loads,  even 
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when  the  resultant  load  on  the  perimeter  of  the  hole  is  the  same.  Fatigue  lifetimes  will 
therefore  be  strongly  dependent  on  the  load  distribution  since  cracks  are  short  for  most 
of  the  lifetime. 

2  BASIC  THEORY 

The  hole  is  defined  in  plane  polar  co-ordinates  (p,6)  by  0  <  p  <  R  and 
0  *»  ('  <  2ir  ard  the  crack  by  R  p  <  a  along  0  «  0  (see  Fig  1);  and  the  crack 
length  i  measured  from  the  edge  of  the  hole  is  given  by  l  •  a  -  R  .  If  the  component 
were  uncracked,  the  applied  loads  would  cause  stresses  o  and  o  .  along  the  crack 

C?D  P“ 

site;  which  can,  in  general,  be  represented  as 

°oe(p,0)  "  P0fi(r)  and  cp9(o,0)  "  Pof2(r)  (l) 


for  R  <  p  <  a  ,  where  o„„  and  a  .  are  the  normal  and  shear  stresses  respectively, 

00  p  o 

r  =  p/R  ,  and  is  a  constant  with  the  dimensions  of  stress.  The  loading  functions 

fj(r)  and  f2(r)  for  the  unerseked  component  are  needed  to  obtain  stress  intensity 
factors  for  the  crack  component. 

2 

It  has  been  shown  that  and  K.^  can  be  expressed  in  terras  of  two  functions 

P,(c)  and  P2(c)  ,  where  c  =  a/R  ,  as  follows 


a 


<C  -  I)  ‘1 


P . ( c )  and 


II 


Si 


(c  -  I)  ‘2 


P,(c)  , 


(2) 


where  ■  p^/iTt  •  The  functions  and  P.,  are  solutions  of  the  singular  integral 

equations. 


Pj(t) 


V  m  -  TTTt  -  0 


and 


P.(D  •  0  , 


-  f.(r),  (I  <  r  <  e)  (3) 


(4) 


for  i  *  1  and  2  ,  where  the  kernel  functions  k.(r,t)  are  given  by, 


Vr.«>  •  -l'-al4  -  - 


and 


t(I  -  rt)’ 


(t  -  t2)2 


(I  -  rt) 

U 


2  (1  -  rt) 


r2t 


t2) 


t(l  -  rt)  tO  -  rt) 


i '  mnn 


(3) 

(6) 


It  should  be  noted  that  the  kernel  function  k j  ( r , t )  used  to  obtain  for  a  loading 

which  is  not  symmetric  about  the  crack  line  is  identical  with  the  kernel  function  k(r,t) 
used  in  Ref  l  for  symmetrical  loading,  because  of  this,  Kj  for  non-syntaetrical  load¬ 
ings  can  be  obtained  from  solutions  for  symmetrical',  loadings,  see  Appendix.  There  is  no 
equivalent  relationship  for  Rj»  ,  since  K,,  is  identically  aero  for  syaaitrieal 
loadings. 


o  BO 
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The  function  Pj  is  related  to  the  displacement  u^  of  the  crack  surfaces,  and 

P  is  related  to  the  displacement  u  : 

»  P 


ufl(p,0)  -  uQ(p,2ir) 


4(1  -  OPqR  r  Pj  (t)dt 


v  ;poK  f 

e  I: 


,  1  <  r  <  c  (?) 


(c  -  t)(t  -  I) 


u  (p  ,0)  -  u  (p,2it)  - 

P  P 


4(1  -  u2)p-R  Cr  P0(t)dt 


;^0  f  2 

?:  I/tr^ Ti 


(t  -  I) 


,  I  <  r  <  c  (8) 


where  v  is  Poisson's  ratio  and  E  is  Young's  modulus. 

3  NWffiRICAl.  PROCEDURE 

The  equations  (3)  and  (4)  are  reduced  to  a  system  of  simultaneous  linear  algebraic 
equations  by  using  a  procedure  developed  by  Erdogan  and  Cupta^.  l.et 


and  v. 

1 


j  ■  1 , 2, . . . ,w  ; 


and  let 


\  m  He  -  \\  ♦  Jit  ♦  l>, 


k  **  I ,ra 


r .  *  J (c  -  l  )v  .  ♦  i  (c  ♦  I ) , 

J  J 


j  •  1,2 . (in  -  1 )  .  HD 


Equations  ())  and  (4)  can  now  be  reduced  to  the  following  linear  algebraic  system: 


m 

»  *  J 


~  Mr.)  . 

*  J 


C<-  - »  • 


with  j  a  I,  2,  ...,  (m  -  1)  and  «  a  1 ,  2  .  The  solutions  P(t  )  of  equation  (12)  are 
used  in  a  Causs-Chebyshev  interpols *io  .  f^-mula  to  determine  the  stress  intensity  factors 


as  follows; 


m  n — *" 

/5  V  kl'* 


r-x 

j  ‘  \  W  *  1  *  **"  ' 


where  tor  convenience  and  K,  «  Kjj  .  The  number  ot  terns  m  is  chosen  large 

enough  to  ensure  that  equation  (13)  converges  and  gives  values  of  K.  to  the  desired 
aeeuraey. 
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4  LOADING  FUNCTIONS 

In  order  to  obtain  the  solution  to  equation  (12),  the  loading  functions  f^(r)  and 
f-,(r)  must  be  known.  These  are  obtained  from  the  stress  fields  along  the  crack  site  in 
the  uncracked  configuration.  In  what  follows,  the  stress  fields  which  are  considered  arc 
due  to  either  uniform  applied  stresses  remote  from  the  crack  or  radial  forces  on  the  hole 
perimeter.  Residual  stress  fields  which  may  also  be  present  in  practical  components  can 
be  included  in  f|(r)  and  f9(r)  .  In  fact,  the  residual  stress  fields  which  often 
occur  around  holei  may  have  a  large  effect  on  the  stress  intensity  factor  particularly  at 
short  crack  lengths. 

4. I  Uniaxial  tensile  stress 

2 

The  stress  intensity  factors  for  this  configuration  have  been  evaluated  by  using 
loading  functions  obtained  from  Sokolnikof f ^ ,  ie  for  a  uniform  uniaxial  tensile  stress 
acting  at  an  angle  if  to  the  crack  line  (see  Tig  la) 


f,(r)  =  ~  i(l  ♦  cos  2$ 

'  r  '  r  ' 

and 

f2(r)  »  i  (l  ♦  -4j- - sin  2<f  . 

r  r 

In  this  case  the  constant  pQ  in  equation  (I)  is  equal  to  the  applied  stress  o 
4.2  Point  load 

The  loading  functions  needed  for  a  radial  force  F  per  unit  thickness  acting  on 
the  hole  perimeter  at  an  angle  4  to  the  crack  line  (see  Tig  lb)  have  been  derived*  from 
the  stress  functions  given  by  Green  and  Zerna  .  The  constant  p  is  given  by 
Py  «■  F / ( 2 R )  .  The  functions  can  be  writteu  as  follows: 

£|(r)  •  U,(*)  *  «g2U)  and  f2(r)  •  g.,(r)  ♦  ag^(r)  (15) 

where  a  =  Y( | (plaue  strain),  or  a  =  — (plane  stress)  .  i  16) 


and 


Sj  (r) 


jJ-L 

»  1 

Ir 


,  4  2r  cos  4  (1  ♦  4r2  -  rS  cos 

♦  -  -  A1  1  ■  ■  ■ 11  • 


rX 


J  * 


(I?) 


e,(r> 


*  r 

~s—  eos  4  , 


*r 


(13) 


(I  -  rV  .  . 

- - -  Sin  4  , 


(1$) 


S4(r) 


I  -  r  * 

- rr-  am  4  , 


*r 


(20) 


f 


o 

% 

ir 


i 


i 
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where 


X^  -  1  -  2r  cos  $  ♦  r^ 


If  we  write 


f’,U)  -  Qj  (t)  ♦  c»Q2(t) 

P2(c)  -  q3(t)  ♦  «q4 ( t )  , 


where  Q^(s  *  l, 2, 3, 4)  is  independent  of  Poisson's  ratio,  and  introduce  the  functions 


Mj  (r,t) 


M3(r,t) 


K2(r,t) 


M4(r,t) 


T^T  ♦  k|(r,t) 


T  *  Vr'°  * 


then  equations  (3)  and  (4)  are  equivalent  to 


f  V0Ms(r'0 


-  8 .Cf)  » 


1  <  r  <  c  (24) 


Q  0)  *  0  (25) 

for  s  »  1,  2,  3,  4  .  Therefore,  by  the  method  of  Erdogan  and  Gupta  (see  section  3),  the 
linear  set  of  equations  to  be  solved  for  Qs(0  is 


iu 

iXvyvyV  B  “  My 


ii« 

5> 


»\<vjri^  '  0 


for  j  «  I,  2, 
given  by 


, ,  (a  -  I)  and  s  *  1,2,  3,  4  .  The  stress  intensity  factors  are 


Kj  »  K(l)  +  uK(4)  and  X.,  »  K(3)  ♦  M4) 


where 


m 


'Hr—W 


for  s  *  I,  2,  3,  4  . 


4 . 3  Are  of  unifora  pressure 

The  stress  function  derived  by  Ureeu  and  ZenM  has  been  used3  to  obtain  the  load¬ 
ing  functions  f  (r)  and  f^(r)  for  the  case  of  a  hole  loaded  by  a  uniform  pressure  on 
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an  arc  of  the  perimeter;  the  arc  on  which  the  pressure  acts  subtends  an  angle  28  at  the 
centre  of  the  hole  and  is  synaetric  about  the  line  0  ■  ir/2  (see  Fig  1c).  As  in  the 


previous  section,  we  write 


fj(r)  “  g]  (r)  ♦  ctg^if)  and  f7(r)  -  g3(r)  ♦  ag^(r)  , 

I  -1/  sin  28  \  P  f  1  V  1  1  \cos  8 

g  (r)  -  -  tan  (— - )  ♦  —  -  [r  ~  7J(—  "  ~ T"  * 

‘  '  r  ♦  cos  28  r  '  ''r,  r/ 


2  1 


g.,(r)  "  0  . 


1/  l\  r  -  sin  B  _  r  ♦  sin  6 

2ff  \  ~  T/  2  ~  2 

L  r,  J 


l/l  sin  0 

«4(r)  =  l(-  '  ‘j  — —  * 


2  2 
♦  2r  sin  8  ♦  r  and  r.; 


2r  sin  8  ♦  r* 


The  stress  intensity  factors  are  calculated 1  by  following  the  steps  of  section  4.2 
f  i  or  equation  (22)  onwards,  but  with  the  following  limitation:  since  g7(r)  =  0  then 
Q  ,  ( t )  “0  and  M.,(r,c)  “  0  and  hence  »  0  . 

4.4  Cosine  distribution  of  pressure 

In  thic  case  the  pressure  p ( 0 )  is  distributed  on  the  perimeter  of  the  hole 


according  to 


p(0)  »  pQ  uos^  ~  8)  =  PQ  8  ,  0  ^  S)  <4  w  , 


and  is  zero  otherwise. 

The  loading  functions  for  any  arbitrary  distribution  of  pressure  can  be  derived,  using 
the  principle  of  superposition,  from  those  for  a  radial  force  acting  on  the  perimeter 

(see  section  4.2).  For  the  case  of  a  cosine  distribution  of  pressure  the  loading 
3  K 

runctions  are  given  by 

f((r)  »  Ij  ♦  and  f7(r)  »  tj  ♦  ,  (3 


Tt 

l  »  1  /  g  sin  $d$  , 

s  J  s 


a  -  1.2. 3, 4  ;  (36) 


g  for  s  •*  1,  2,  3,  4  are  the  loading  functions  lor  a  point  load,  given  by  equations 
(17)  and  (20).  The  integrals  I  are  evaluated,  for  s  »  I,  2,  3,  4  ,  in  Appendix  8  of 
Ref  1;  the  expressions  for  fj(r)  and  f,(r)  become 


i(,i  •  ^[5  *  7  *  (ri  *  7  -  tMt^-!)] 


»«0 


990 
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and 


Mr) 


Y(»  -  r  ) 
4r3 


(38) 


where  y  -  •?  -  1  . 

W:.th  the  following  definitions, 


and 


h((r)  »  f  j (r)  and  h  2 ( r )  =  -f2(r) 


Q j (r)  -  P j (r)  and  Q2(r>  ”  -  r,(r) 


equations  (3)  and  (4)  become 


(39) 


Qj(t) 


✓  fc  -_t)(T-  I ) 


— ♦  k.(r,t)|  -  -  h.  ( t)  ,  1  <  r  <  c  , 


(40) 


with  Q^(l)  -  0  ,  for  i  ■  1 ,  2  .  The  stress  intensity  factors  are  given  by 


/2 


(c  -  13 


T  I 

Q.(c)  and  — — 


f/2 


(c  -  1)  ^2 


0>(c)  . 


(41) 


The  numerical  procedure  to  solve  equation  (40)  and  hence  obtain  and  is  as  out¬ 

lined  in  section  3. 

5  RESULTS 

Solutions  of  the  singular  integral  equations  (3)  and  (4)  in  section  2  have  been 
obtained  and  hence  the  stress  intensity  factors  evaluated,  by  using  the  numerical  pro¬ 
cedure  outlined  in  section  3  and  the  loading  functions  developed  in  section  4. 

5. I  Uniaxial  tensile  stress 

Values  of  the  opening-mode  stress  intensity  factor  K^/K^  are  Riven  in  Table  I  as 
a  function  of  the  crack  length  i/R  for  various  values  of  the  angle  *  between  the 
direction  of  the  stress  and  the  eraek  line.  The  results  for  $  of  0°  and  90°  agree  with 
those  for  biaxial  and  uniaxial  stress  given  in  Ref  I.  Plots  of  K  /KQ  against  t/R  are 
shown  in  Pig  2  for  0°  *=  i  <  90°  .  For  other  values  of  $  the  stress  intensity  factor 
KjU)  is  obtained  from  the  following  syeaaetry  relations: 

R,(-  i)  *  -  $)  KjM  ,  0°  <  i  <  180°  .  142) 

Negative  values  of  Rj  inply  eraek  elosure  and  are  therefore  not  generally  valid  but 
they  ean  be  superposed  with  values  for  other  stress  fields  provided  that  the  resultant 
stress  field  causes  the  eraek  to  remain  open  along  its  whole  length. 


T 


Values  of  the  sliding-node  stress  intensity  factor  are  8*,von  *n  2 

and  plotted  in  Fig  2.  The  sytanetry  relations  arc 


K  (-  *)  -  K  (180°  -  $)  -  -  K  (♦) 


K1t(90°  -  «>  *  Kn(0)  ,  0°  <  <j>  «  180 


.  } 


5 . 2  Point  load 

The  components  K^/K^  and  of  the  opening-mode  stress  intensity  factor 

Kj/'Kq  are  given  in  Tables  3  and  4,  respectively.  Corves  of  these  two  components  are 
plotted  as  a  function  of  i/R  in  Fig  3  for  various  values  of  <f>  .  For  clarity  of  pre¬ 
sentation  is  plotted  for  0°  4.  f  <S  180°  and  K^/K^  for  90°  *»  $  <  180°  . 

Results  tor  other  values  of  $  can  be  obtained  from  the  following  symmetry  relations: 


K(1)(-4>  =  K(l)(» 


K.(~\-  *)  -  K(i)(.j)  -  -  K(2)(180°  -  *)  ,  0°  <  * 


0  a  <.  iso0  .  1 


The  components  K^/Kq  and  K^4*/Kq  of  the  sliding-mode  stress  intensity  factor 
Kn/R0  are  given  in  Tables  5  and  b  respectively.  These  are  plotted  as  a  function  of 
i/ R  in  Fig  4  for  various  values  of  $  .  The  symmetry  relationships  are 


K*  ,}(-  *) 


~  K(J)(+) 


K(4,(180°  -  *)  =  K(4)U)  =  -  R(4>(-  *>  .  0°^*<180°. 


In  general  and  K.^  are  functions  of  Poisson’s  ratio  v  ,  for  this  loading, 

since  from  equation  (27)  we  obtain,  for  plane  straiu. 


<  3  -  4v)  K  ^  ~ } 

2Ti"-V>  K„ 


(3  -  4v) 

TO  -  v)  *>„ 


Values  of  and  have  been  calculated  for  v  -  0.1  and  are  shown  in  Figs  5 

and  a  respectively  as  a  function  of  t/%  for  various  values  of  the  angle  f  .  For 

o  t  j )  ,  ... 

}  —  VO  the  eotapuneat  R  —  0  and  the  opcaing'taodc  a[*'ca4  intensity  factor  is 

independent  of  Poisson's  ratio,  and  is  otto  half  the  vale-  retained  previously*  for  two 

equal  and  opposite  forces  acting  on  the  hole  perimeter. 

-  CJ 

For  i  **  11  crack  closure  can  oeeur  for  some  crash  lengths,  and  if  it  does  occur 
the  values  of  R^  calculated  by  the  present  methods  are  invalid;  if  closure  occurs  aC 
the  tip  Rj  will  be  negative.  The  values  of  obtained  can  be  used  only  if  another 


**!£> 
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stress  field  is  superimposed  on  that  due  to  the  point  load,  such  that  the  crack  faces  do 
not  touch  at  any  point  along  the  length  of  the  crack. 

5 -  3  Arc  of  uniform  pressure 

Values  of  the  opening-mode  stress  intensity  factor  K^/K^,  are  given  in  Tables  7aib 
and  are  plotted  in  Fig  7  as  a  function  of  the  angle  ?>  for  various  values  of  the.  crack 
length.  The  normalising  constant  Ky  is  given  by 


Kp  "  kq  *in  B  -  pQ-' 


it l  sin  8 


It  is  convenient  to  use  Ky  since  it  is  proportional  to  the  resultant  for'’  K  sin  ?• 

acting  on  the  hole  perpendicular  to  the  crack-line.  1’or  6  •  0  the  re1-  .'.s  agre--  with 
K^^/Kq  for  a  point  force  acting  at  an  angle  |  »  90°  to  the  crack*  ’  '  ,ie  (Table  3);  for 
ii  “  90°  the  values  of  /K^_,  are  one  half  of  the  values  obtained  pra.-.ously  for  a  uni¬ 
form  pressure  acting  over  the  whole  perimeter. 

The  values  of  K./K.,  given  in  TabLes  7aAb  and  Fig  7  are  tor  C°  9  ft  w  90°  ;  values 

l  1' 

for  d  >  90  can  be  obtained  as  follows: 

H  (B)  ♦  K-l ( 1  80°  *  B )  =  21^(90°)  ,  0°  *  ft  <  130°  .  (48) 

The  two  components  K^/Ky  and  K^/Ky  of  the  sliding-mode  stress  intensity 
•  ( 4 ) 

factor  are  tabulated  in  Tables  SaAb;  the  function  K  is  a  function  of  crack 

length  only,  independent  of  r.he  angle  8  .  Values  of  Vk.  are  plotted  as  a  function 

o  *  (  '  ) 

of  B  ,  for  0  *«  B  **  90°  ,  in  Fig  8a  for  various  crack  lengths.  The  component  K  "*  /KB 

r 

is  plotted  in  Fig  8b  as  a  function  of  t/U  . 

The  symmetry  relationship  for  is  as  to l  lows: 


KU)(B)  *  K(J)U80°  -  B)  , 


0°  *  g  'S*  ISO1*  .  149) 


4,4  Cosine  distribution  of  pressure 

Values  of  t£j/ky  and  f  ^(K^/Ky)  ate  tabulated  in  Tabic  9  for  the  cosine  dietri* 
button  of  pressure;  for  this  ease  Sh,  is  given  bv 


■?  *  31  ’*£ 


Where  F  is  the  resultant  forec  acting  perpendicular  to  the  crack-line  and  is  given  by 


thus 
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The  opening-mode  stress  intensity  factor  is  plotted  as  K^/K^,  against  2./R  in 
Fig  9.  Values  of  Kj/Kp  f°r  two  other  lead  distributions  having  the  same  resultant 
force  perpendicular  to  the  crack-line  are  also  plotted  in  Fig  9;  the  two  distributions 
are  a  point  force  acting  at  an  angle  <f>  *»  90°  (see  Fig  5)  and  a  uniform  pressure  acting 
on  a  semi-circular  arc  of  the  perimeter  (see  Fig  7).  Comparison  of  these  three  curves 
shows  that  although  the  resultant  force  is  the  same,  the  different  load  distributions 
lead  to  differences  in  the  stress  intensity  factors,  and  these  differences  are  greatest 
when  the  cracks  are  short. 

Fig  9  also  shows  a  plot  of  y  ^K^/Kp  against  i/R  . 

6  DISCUSSION 

The  method  described  in  this  Report  has  two  main  advantages;  it  is  efficient  in 
computing  requirements  and  it  is  accurate  at  short  crack  lengths.  The  largest  number  of 
simultaneous  linear  equations  needed  to  obtain  results  accurate  to  ~0.1”  was  24.  Hsu^ 
who  has  solved  the  crack  problem  for  uniaxial  tensile  stress  only  (see  Fig  la)  required 
between  two  and  six  times  as  many  equations.  His  method  which  uses  conformal  mapping 
techniques  becomes  less  accurate  for  short  cracks,  in  contrast  to  the  present  method  for 
which  the  stress  intensity  factor  tends  towards  the  correct  limit,  for  a  crack  at  the 
edge  of  a  hole,  as  5./R  tends  to  zero.  This  ensures  that  values  of  stress  intensity 
factors  are  accurate  for  short  cracks,  and  this  is  a  necessary  requirement  for  the 
estimation  of  fatigue  lifetimes  to  be  reliable. 

This  limiting  value  can  be  expressed  in  terms  of  the  loading  functions, 

Lt  )K  }  -  1.122c  (1, 0)/n  -  1.122p  f  (1)^1  (53) 

1/  R+0 

where  the  factor  1.122  is  the  usual  edge  correction  factor.  It  therefore  follows  from 
the  definition  of  Kq(=  Pq*^)  that 

Lt  fej  -  !.122f  (!)  .  (54) 

i/R-*0(K0) 

For  a  point  force  (<p  =  90°),  f  ((1)  =  2/ir  from  equations  (15),  (17)  and  (21);  for  an  arc 
of  pressure  (B  =  90°),  f  f  ( 1 )  =  1/2  from  equations  (29)  and  (30);  and  for  a  cosine  dis¬ 
tribution  of  pressure,  fj(l)  =  2/u  from  equation  (37). 

Comparison  of  the  stress  intensity  factors  for  the  point  force,  the  arc  of  pressure 
and  cosine  distribution  of  loads  snows  (see  Fig  9)  that  the  limiting  values  of  K^/K^, 
are  1.122  x  (2/n),  l.i22  x  (1/2)  and  1.122  x  (8/tt^)  respectively;  they  are  in  the 
ratios  1  :  0.78  :  1.27  .  Thus  even  when  the  three  load  distributions  have  the  same 
resultant  force  F  ,  the  values  of  for  short  cracks  differ  significantly  from  each 

other.  Because  of  this  the  precise  load  distribution  must  be  known  if  practical 
estimates  of  fatigue  life  and  residual  strength  are  to  be  accurate.  This  is  particularly 
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important  in  the  calculation  of  fatigue  life  because  the  growth-rate  of  cracks  in 
fatigue  is  proportional  to  Kn  where  n  is  of  order  2  to  4. 

It  is  most  important  that  the  load  distribution  near  the  root  of  the  crack  is 
known  precisely,  since  small  changes  in  load  in  this  region  result  in  large  changes  in 
K  .  This  can  te  seen  from  Figs  3  and  6  where  both  and  change  rapidly  as  the 

angle  $  between  the  force  and  the  crack  becomes  small.  A  similar  effect  is  seen  in 
Figs  7  and  8,  where  rapid  changes  occur  at  small  crack  lengths,  as  the  arc  over  which  the 
uniform  pressure  acts  approaches  the  crack  (B  «  90°). 

The  stress  intensity  factors  and  ,  obtained  in  this  Report  are  for  a 

single  crack  at  the  edge  of  a  hole  with  an  arbitrary  distribution  of  load  on  the  hole 
perimeter.  Values  of  the  only,  for  two  cracks  of  unequal  length  at  •>  »  0  and  it  , 

can  be  obtained  for  arbitrary  loads  by  using  the  principle  of  superposition  (see 
Appendix)  and  the  method  developed  in  Ref  1.  A  second,  more  approximate  method',  can 
also  be  used  for  obtaining  for  two  cracks  from  that  for  one  crack  provided  that  the 

load  distributions  on  the  perimeter  is  symmetric  about  9  *  tt / 2  .  Although  only  the 
opening-mode  stress  intensity  factor  can  be  obtained  from  both  procedures  this  is  often 
adequate  for  practical  structures,  cince  is  usually  much  more  than  and  the 

fracture  behaviour  is  dominated  by  cho  opening-mode. 

The  method  used  in  this  Report  to  evaluate  stress  intensity  factors  for  the  four 
•  ■ »  hi  nlar  loading  distributions  shown  in  Fig  I,  can  be  extended  to  any  arbitrary’  loading 
on  U»u  hole  perimeter  by  a  technique  using  Green's  functions7.  This  technique  was  usee 
in  section  4.4  to  obtain  the  loading  functions  f  (r)  and  f ..,  ( r  )  for  the  cosine  distri¬ 
bution  of  pressure.  For  an  arbitrary  pressure  distribution  p ( t* )  on  the  hole  perimeter, 
which  can  be  written  as 

p(h)  =  pQh(0)  ,  (35) 

equations  (35)  and  (36)  become 

f  (r)  =  I ,  ♦  al  ,  and  f  ..(r)  »  1..  ♦  nl  (56) 

1  i  i  2  3  4 

2n 

where  “  j  J  gg U)h ($)d$  ;  (57) 

0 

the  Green's  functions  g  ($)  *  for  s  “  1.2,  3  and  4  ,  are  the  loading  functions  for  a 
point  load,  and  are  given  by  equations  ('7)  to  (20).  After  evaluation  of  I  ,  analytic¬ 
ally  or  numerically,  the  procedure  of  section  3  can  be  followed  to  obtain  the  stress 
intensity  factors. 

The  stress  intensity  factors,  obtained  in  section  5.2  for  a  point  load,  may  be 
used  directly  as  Green's  functions7.  Let  che  force  acting  on  the  perimeter  between 
6  »  i]>  and  6  “  $  +  d<}>  due  to  the  pressure  p(6)  be  dl'($)  ;  it  is  given  by 

dF(<&)  “  pti>)Rd$  .  (58) 


u 


Lei:  the  stress  intensity  factor  (I  or  II)  due  to  dF($)  be  dK($)  ;  it  is  given  by 


dK($)  ■  K($)dF($)  , 


(59) 


where  KU ( )  is  the  factor  for  a  unit  force  acting  at  9  -  #  .  Since  in  general  the 
stress  intensity  factor  is  proportional  to  the  applied  force  F  it  follows  that 


KU(*) 


(60) 


where  K  (■{')  is  the  factor  for  a  force  F  acting  on  the  hole  perimeter  at  0 
Equation  (54)  can  be  written 


KUU) 


K*U) 


■U  it) 

2R  F  r-r 

nr1 


>WkF(/,)\ 
2R  V  Kn  ; 


(61) 


where  (K  ($)/KQ)  can  be  either  K^/K^  or  K  /K^  obtained  in  section  5.2. 

The  principle  of  superposition  states  that  the  stress  intensity  factor  k'*  duo  to 
p(t)  for  all  is  given  by  summing  all  the  dK(o)’s,  %t 


/ 


dK(  i£) 


(62) 


Substitution  of  equations  (55)  and  (58)  to  (61)  into  equation  (62)  gives 


K 


P 


(63) 


Thus  (K^  (^)/K(j)  is  the  Green’s  function  for  the  stress  intensity  factor  for  an  arbitrary 

distribution  of  pressure  and  can  be  obtained  from  the  results  for  a  single  point  force. 

The  integral  in  equation  (63)  will,  in  general,  have  to  be  evaluated  numerically.  If  the 

loading  distribution  is  determined  exper tiacntdl ly  so  that  h($)  is  known  only  at  discrete 

points,  then  the  integral  in  equation  (63)  would  be  replaced  by  aa  appropriate  summation. 

in  order  to  ensure  accurate  results,  particular  care  muiet  be  taken  when  3  is  near  seru 
F 

as  K  ($)  varies  rapidly  (see  Figs  5  and  6). 

7  COhCLUSlOMS 

(l)  A  transform  method  has  been  developed  for  obtained  opening-mode  and  sliding-mode 
stress  intensity  factors  for  a  crack  at  the  edge  -->i  a  circular  hole;  the  method  is  both 
accurate  and  efficient  in  computer  usage.  The  method  was  developed  so  that 

(a)  the  stress  intensity  factor  for  an  arbitrary  distribution  of  pressure  in  the 
hole  can  be  obtained  from  results  for  a  radial  point  force; 


vso 


(b)  both  externally  applied  loads  and  residual  stress  fields  can  be  considered; 

(c)  it  can  be  extended  to  calculate  opening-mode  stress  intensity  factors  for  two 
cracks . 

(2)  Application  of  the  method  to  a  short  crack  at  the  edge  of  a  hole,  showed  that  the 
stress  intensity  factor  varies  significantly  with  small  changes  of  load  distribution  on 
the  hole,  even  if  the  total  load  on  the  hole  is  the  same;  this  implies  that  estimates  of 
fatigue  lifetimes  may  be  seriously  in  error,  unless  the  load  distribution  in  the 
structural  component  is  known  accurately. 
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Appendix 

RELATIONSHIP  BETWEEN  STRESS  INTENSITY  FACTORS  FOR  SYMMETRIC 
AND  NON-SYMMETRIC  LOAD  DISTRIBUTIONS 
(see  section  2) 

By  using  the  principle  of  superposition,  any  load  distribution  which  is  symmetric 
with  respect  to  the  crack  line  can  be  represented  as  the  sivm  of  two  asymmetric  loadings, 

In  Fig  10,  configuration  A  represents  a  hole  with  a  crack  at  its  edge  along  the  0=0 

A 

axis;  the  hole  is  loaded  symmetrically  about  the  crack  line,  so  that  the  pressure  p 
on  the  perimeter  is  given  by 


P  =  P (0)  » 


P("  8)  , 


0  <  0  ^  7T 


0  >  0  '-2  -IT  . 


Configuration  B  represents  a  similar  hole  and  crack  with  loading  for  0  <  0  ^  it  only,  so 

g 

that  the  pressure  p  on  the  perimeter  is  given  by 


P  =  p(9)  , 


P  =  0  , 


0  «  6  it 


o  >  e  >  -it  . 


Configuration  C  also  represents  a  similar  hole  and  crack,  but  is  loaded  for  0  >  0  >  -n 

C 

only,  so  that  the  pressure  p  o.i  the  perimeter  is  given  by 


P  =  0  , 


Pc  =  p(-  o)  , 


0  <  0  <  n 


0  **  0  >  -IT 


The  stress  intensity  factors  for  these  three  configurations  can  be  written  as 


follows: 


for  A 


Kt  -  K*  and  Kn  =  ; 


for  B 


K!  ‘  K?  and  KII  *  KII  > 


for  C 


Kj  =  and  Kn  =  ; 


From  the  principle  of  superposition  it  follows  that 


il  K11 


From  symmetry  considerations  it  follows  that 


kj  «  hj  and 
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hence 


and 


=  0 


(A-9) 


Thus  it  follows  from  equation  (A-9)  that  the  opening-mode  stress  intensity  factor 
for  an  asymmetric  distribution  of  loads  can  be  obtained  from  K.^  for  a  symmetric  distri¬ 
bution;  but,  this  is  not  so  for  K  since  it  is  identically  zero  for  symmetric  loadings. 

Because  this  argument  applies  to  configurations  with  one  crack  (along  9  =  0)  or  two 
cracks  (along  0=0  and  0  =  it),  the  method  developed  previously',  for  obtaining  for 

two  symmetrically  loaded  cracks  of  unequal  length,  can  be  used  to  obtain  for 

asymmetrical  loadings.  The  factor  cannot  be  determined  in  this  way;  but  for  many 

practical  loadings  (K^  |  ^  ,  so  that  the  fracture  of  the  cracked  component  is 

dominated  by  the  opening-mode  stress  intensity  factor. 


036 


Table  1 


OPENING-MODE  STRESS  INTENSITY  FACTOR  K_/K„  FOR  A  CRACK 


AT  THE  EDGE  OF  A  HOLE:  UNIAXIAL  TENSILE  STRESS 


30°  45°  60°  75°  90° 


\  * 

0° 

n/ik 

0.01 
0.02 
G.04 
0.07 
0. 10 
0.15 
0.20 
0.30 
0.40 
0.50 
0.60 
0.80 
1.00 
1.50 
2.00 
3.00 
4.00 
5.00 
6.00 
8.00 
10.00 
20.00 


-1.080 
-1 .040 
-0.967 
-0.869 
-0.783 
-0.664 
-0.567  I 
-0.421 
-0.320 
-0.247 
-0.194 
-0.123 
-0.080 
-0.030 
-0.010 
0.002 
0.005 
0.005 
0.005 
0.004 
0.003 
0.G02 


0.013 

0.026 

0.049 


-0.545 

-0.448 

0.106 

0.141 

-0.370 

0.169 

-0.253 

0.207 

-0,172 

0.232 

-0. 115 

0.247 

-0.073 

0.256 

-0.019 

0.265 

0,013 

0.266 

0.048 

0.260 

0.059 

0.250 

0.064 

0.234 

0.063 

0.223 

0.061 

0.215 

0,060 

0.210 

0.057 

0.202 

0.055 

0.  197 

0.051 

0.187 

1 . 106 
1 .092 
1.065 
1 .028 
0.994 
0.946 
0.904 
0.836 
0.783 
0.740 
0.706 
0.652 
0.613 
0.549 
0.310 
0 , 466 
0.441 
0.425 
0.414 
0,400 
0.391 
0.373 


2. 199 
2.158 
2.081 
1 .976 
1 .883 
1.750 
1.639 
1 .464 
1.334 
1.234 
1.155 
1.040 
0.959 
0.838 
0.770 
0.698 
0.b59 
0.635 
0.619 
0.598 
0.585 
0.559 


2.999 
2.938 
2.824 
2.670 
2.534 
2.339 
2.177 
i  .924 
1.737 
1 .595 
1.484 
1.323 
1  .213 
1 .050 
0.96 1 
0.868 
0.819 
0.789 
0.769 
0.743 
0.727 
0.695 


3.292 
3.224 
3.096 
2.924 
2.772 
2.555 
2.374 
2.092 
1.885 
1.728 
1 .605 
1.427 
1.306 
1.127 
1.031 
0.930 
0.878 
0.846 
0.824 
0.796 
0.779 
0.744 


E 


Table  2 

SLIDING-MODE  STRESS  INTENSITY 


FACTOR  K„/K„  FOR  A  CRACK 


AT  THE  EDGE  01  A  HOLE 


UNIAXIAL  TENSILE  STRESS 


\  * 

0° 

15° 

0.01 

0.000 

0.013 

0.02 

0.000 

0.026 

0.04 

0.000 

0.050 

0.07 

0.000 

0.081 

0.  10 

0.000 

0.  109 

0.15 

0.000 

0.147 

0.20 

0.000 

0.178 

0.30 

0.000 

0.223 

0.40 

0.000 

0.252 

0.50 

0.000 

0.271 

0.60 

0.000 

0.284 

0.80 

0.000 

0.297 

1 .00 

0.000 

0 .  300 

1 . 5C 

0.000 

2“  t 

2.00 

0 . 000 

0.283 

3.00 

0  .  JOO 

0. 2o 3 

4.00 

0.000 

0.248 

5.00 

0 . 000 

0.237 

6.00 

0.000 

0.229 

8.00 

0.000 

0.216 

10.00 

0.000 

0.204 

20.00 

0.000 

0.144 

0.023 
0.045 
0.086 
0.141 
0.189 
0.255 
0.306 
0.386 
0.437 
0.470 
0.492 
0.514 
0.520 
0.51 0 
0.490 
0.455 
0.429 
0.4  1  1 
0.39? 
0.  374 
0.  362 
0.3  36 


0.02  7 
0.052 

n.  ioo 

0. 163 
0.218 
0.295 
0.356 
0.446 
0.504 
0.543 
0.568 
0.593 
0.600 
0.384 
0.566 
0.525 
0.496 
0.474 
0.458 
0.432 
0.418 
0.  388 


Table  .3 

COMPONENT  OK  THF.  OPENING-MODE  STRESS  INTENSITY  FACTOR  FOR  A 


CRACK  AT  THE  EDGE  OF  A  HOLE:  POINT  LOAD 


15° 

: 

30° 

\ 

,  -  O 

i  > 

0.579 

0.672 

0.690 

0.4  AH 

O.t.  53 

0.666 

t).  32  ! 

0 . 49  5 

0.644 

0.705 

0.555 

0.621 

-0.009 

0.483 

0.579 

-0.190 

0.417 

0.540 

-0,039 

0.355 

i 

0.503 

-0.583 

1 

0.222 

0.420 

-0.688 

0.118 

0.  549 

-0.688 

-0.021 

0.238 

-0.601 

-0.096 

0.  IbO 

-0.,  1? 

-0. 148 

0.065 

-0.  5  u 

-0.1 40 

0.01  / 

-0.765 

-0.140 

-0.008 

-0. 1 79 

-0 .  1  1  1 

-0.031 

-0.1  58 

-0.091 

-0.053 

-0.097 

-0.068 

-0.054 

-0.0/6 

-0.036 

-0.051 

-0.054 

-0.041 

-0.023 

-0.043 

-0.055 

| 

-0.031 

1 

C 

c 

-0.028 

-0.013 

-0.0P9 

-o.oib 

1 

1 

-0.0 1 1 

u.no  3 
0.000 


0 

l  20 

180° 

0.701 

0.702 

0.690 

0.691 

0.678 

0.680 

0.66  7 

0.670 

0.646 

0.650 

0.626 

0.631 

0.607 

0.614 

0.565 

0.574 

0.529 

0.540 

0.468 

0.484 

0.419 

0.438 

0.34  7 

0.371 

0.295 

0.72  2 

0.255 

0.285 

0.190 

0.22  1 

0.150 

0. 181 

0.104 

0.133 

C 

c 

0.103 

c 

c 

0.074 

0.0  39 

0.037 

0.051 

0 .046 

0.013 

! 

0.02  3 

Tab l c  5 


Table  6 


K0  COMPONENT  OF  THE  SLIDING-MODE  S TRICSS  INTENSITY 


FACTOR  FOR  A  CRACK  AT  THE  EDGE  OF  A  HOLE:  POINT  LOAD 


\  * 
i/iK 

15° 

30° 

45° 

O.OI 

-0.00! 

-0.002 

-0.003 

0.02 
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-0.015 
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1.30 

-0.027 
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-0.033 
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LIST  OF  SYMBOLS 
half  of  total  crack  length 
=  a/R 

Young's  modulus 

for  per  unit  thickness 

force  acting  between  $  and  4«  +  d4> 

loading  functions  (i  =  1,2) 

components  of  loading  functions  (i  =  1,2,3, 4) 

functions  defined  by  equation  (39) 

pressure  distribution  on  the  perimeter  of  a  hole 

integrals  defined  by  equation  (36) 

functions  defined  by  equations  (7)  and  (8) 

constant  pressure  (or  stress) 

summation  variables 

functions  defined  by  equations  (5)  and  (6) 

stress  intensity  factor  SIF  (=  p^/itT) 

opening-mode  SIF 

sliding-mode  SIF 

Kj  and  respectively 

components  of  SIF  (i  =  1,2, 3, 4)  -  see  equation  (27) 
SIF  due  to  resultant  force 
SIF  due  to  dl’ ( d> ) 

SIF  due  to  unit  force  at  ij, 

SIF  due  to  force  F  at  <J> 

SIF  due  to  arbitrary  distribution  of  pressure 
=  a  -  R 

number  of  simultaneous  equations 
functions  defined  by  equation  (23),  (i  =  1,2, 3, 4) 
functions  defined  by  equation  (22),  (i  =  1,2, 3, 4) 
radius  of  circular  hole 
“  p/R 

discrete  value  of  r 


LIST  OF  SYMBOLS  (concluded) 


u 

P 

X 

a 

e 

Y 

e 


v 

P 

o 


♦ 


functions  defined  by  equation  (33) 

variable  of  integration 

discrete  value  of  t 

angular  displacement 

radial  displacement 

function  defined  by  equation  (21) 

-  (3  -  4v)/(  2(1  -  v)]  or  (3  -  v)/2 
angle  (section  4.3) 
a  -  1 

polar  co-ordinate 
Poisson's  ratio 
polar  co-ordinate 
applied  stress 
normal  stress 
shear  stress 

angle  between  force  and  crack  line 
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